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ABSTRACT 

We present a multi- wavelength study of the region towards the GM 2-4 nebula and the 
nearby source IRAS 05373+2340. Our near-infrared H2 1-0 S(l) line observations reveal var- 
ious shock-excited features which are part of several bipolar outflows. We identify candidates 
for the driving sources of the outflows from a comparison of the multi-waveband archival 
data-sets and SED modelling. The SED spectral slope (a^^^^) for all the protostars in the 
field was then compared with the visual extinction map. This comparison suggests that star 
formation is progressing from NE to SW across this region. 

Key words: circumstellar matter - infrared: stars - ISM: jets and outflows - ISM: individual: 
GM 2-4. 



1 INTRODUCTION 

Since 1998, the Byurakan observatory in Armenia has been en- 
gaged in surveys of nebular objects in dark molecular clouds and 
star formation regions. With its 2.6m telescope, imaging in narrow 
band filters and simultaneous searches of Ha emission stars are 
being performed. New HH outflows and star forming regions have 
been found (see Movsessian et al. 2008, and references therein). An 
exploration of the region with GM 1-64 (RNO 53) (Gyulbudaghian 
& Maghakian 1977; Cohen 1980) and the GM2-4 nebulae (Giul- 
budagian & Magakian 1977) at optical wavelengths has lead to the 
discovery of several HH objects as well as numerous Ha emission 
stars (Nikogossian et al. 2009). 

GM2-4 and the nearby IRAS 05373+2349 are located in the 
centre of the [KOY98]183.7-03.6 molecular cloud which has been 
detected in ^^CO by Kawamura et al. (1998). The authors provide 
a distance estimate of ~2kpc. Dobashi et al. (2005) identified a 
peak of optical extinction (TGUH1314 in Fig. 1) from DSSl plates 
which is located about 1.5-2.0' westwards of IRAS 05373+2349 (cf. 
Fig. la). This peak appears to be associated with the Shajn 147 su- 
pernova remnant (Casoli et al. 1986). The distance to the SNR pro- 
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genitor PSR J0538+2817 is obtained from a precise parallax mea- 
surement and is given as 1.3kpc (Chatterjee et al. 2009). These 
distances are smaller than the early estimates of 1.7 kpc based on 
the kinetic velocities of radio lines (eg. Rodriguez et al. 1980), but 
they agree with subsequent radio observations presented by Moli- 
nari et al. (1996) who infer a distance of 1.17 kpc from an updated 
Galactic rotational model (Brand & Blitz 1993) 

GM2-4 was discovered on DSSl plates (Giulbudagian & 
Magakian 1977) and is classified as a Herbig-Haro-Like object 
(HHL 26) by Gyulbudaghian et al. (1978) and later renamed GGD 4 
by Gyulbudaghian et al. (1987). Subsequent near-infrared observa- 
tions (Carballo et al. 1988) were able to separate this object into 
two separate red nebulous objects namely GGD 4a and GGD 4b 
coinciding with the position of GM2-4 (Fig. la). A recent contin- 
uum and narrow-band (Ha and [SII]) optical study by Nikogossian 
et al. (2009) concludes that these objects (GM 2-4 and GGD 4a) are 
reflection nebulae and do not exhibit any Herbig-Haro characteris- 
tics. The near-infrared H2 1-0 S(l) line survey of Simon & Joyce 
(1983) did not find molecular hydrogen line emission, but this is 
likely due to the limited sensitivity of the infrared bolometer used 
in that study. The detection of an H2O maser by Rodriguez et al. 
(1980) and the detection of a cold IRAS point source 05373+2349^ 
about 1.5^ south- west from the centre of GM2-4 (cf. Fig. la) indi- 
cated that active star-formation occurs in the region. 

Casoli et al. (1986) detected high- velocity blue and red- shifted 
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Figure 1. General field including GM 1-64 and GM2-4 centred on IRAS 05373+2349. Panel (a) is DSS2 red image of the region overlaid with a visual 
extinction contours obtained from Dobashi et al. (2005) peaking at TGU HI 3 14. The solid-line ellipse indicates the positional uncertainty of IRAS 05373-1-2349 
and the green asterisks mark the positions of masers. Panel (b) is the 2MASS false colour composite view of the region constructed from J (blue), H (green) 
and Kg (red) bands. Overlaid are the IRAS lOO/^m contours at 50, 70, 90 and 110 MJy/sterad. Panel (c) is Spitzer IRAC colour composite image constructed 
from 3.6/im (blue), 4.5fim (green) and S.O/^m (red) bands. Stars marked with big and small circles are sources reported by Gutermuth et al. (2009) and 
classified as Class I protostars and Class II pre-main-sequence stars respectively. Panel (d) is Spitzer MIPS 24/im band view of the region where the same 
Class I sources are marked with the nomenclature used in our study (see Sect 3.1). 



wmgs m the ^^CO emission which pointed to a bi-polar outflow 
driven by IRAS 05373+2349. Carballo et al. (1988) obtained pho- 
tometric measurements in near- to mid-infrared bands of a red ob- 
ject coincident with the position of the IRAS 05373+2349 using an 
infrared bolometer. The same red object was resolved into an in- 
frared stellar cluster by Hodapp (1994) (cf. Fig. lb). This infrared 
cluster is also associated with 450/im continuum emission (Jen- 
ness et al. 1995) and with an ammonia core (Molinari et al. 1996). 
More recently, 3.4cm and 3.6cm continuum emission and HCO^ 
(1-0), H^^CO+ (1-0) and H^^CN (1-0) line emission were detected 
towards the cluster (Molinari et al. 2002). A fit to the spectral en- 
ergy distribution (SED) of IRAS 05373+2349 from 8 - 1200 /xm by 
Molinari et al. (2008) was interpreted in terms of a zero-age main- 
sequence B5 star of 26 Mq for the entire cluster of stars. The study 
by Gutermuth et al. (2009) using Spitzer addressed the evolution- 



ary stages of individual members in the IR cluster, and found that 
its members were a mixture of Class I protostars and Class II pre- 
main-sequence stars with optically thick disks (cf. Fig.lc). Varricatt 
et al. (2010) reported some very faint near-infrared H2 1-0 S(l) line 
knots around the position of IRAS 05373+2349. 

These studies demonstrate that the area around GM 2-4 is an 
active and on-going star formation region. We have thus decided 
to combine previously published data on the region with our wide- 
field near-infrared survey of the region in the outflow tracing H2 
1-0 S(l) line in order to reveal the extent of outflow activity, iden- 
tify their driving sources and examine the star formation processes 
of the studied field. In Section 2 we present the observations and the 
data analysis. In Section 3 we present the results from our observa- 
tions, discuss the possible outflows and present SED modelling on 



GM2-4 region 3 




OS'^JfrSE" 30' 24' 20* 05' 40"" 32' 2S' 24= 20' 05" 40" 32' 28' 24-' 20' 

Righl Ascension Right Ascension RpgW Ascensmn 



Figure 2. Close-up view of GM2-4 and IRAS 05373+2349 including the region where new molecular hydrogen emission line objects (MHOs) are found. 
Class I sources reported by Gutermuth et al. (2009) are marked with a crosses on all panels. Panel (a) is optical [SII] image where Herbig-Haro objects 
discovered by (Nikogossian et al. 2009) are marked along with the GM2-4, GGD4a and IRAS 05373+2349 error ellipse. Panel (b) is near-infrared H2 1-0 
S(l) line+continuum image where the MHOs are marked with squares. Panel (c) shows Spitzer IRAC colour-composite image constructed from 3.6/im (blue), 
4.5/im (green) and 8.0/im (red) bands with positions of newly detected MHOs marked with squares for a comparison. For this particular case the brightest 
(mostly saturated) stars in the field were truncated in order to stretch the colour scheme and reveal relatively fainter objects. 



suspected driving sources. Finally, in Section 4 we summarise our 
findings and present some future prospectives on the project. 



2 OBSERVATIONS AND DATA ANALYSIS 

The following discussion makes frequent use of the optical obser- 
vations by Nikogossian et al. (2009) - see that paper for a compre- 
hensive discussion of the optical dataset. 



MOPEX software (Makovoz & Khan 2005). Photometric measure- 
ments were done using the aperture photometry package PHO- 
TOM in GAIA^ . Instrumental values (counts/sec) where calibrated 
using the guidelines presented in SSC home page (i.e. appropri- 
ate aperture and pixel-scale correction) and zero-point fluxes by 
Reach et al. (2005) based on Vega-standard magnitudes for 1 DN/s. 
Our results are in excellent agreement with the values presented in 
Gutermuth et al. (2009) differing only by a few percent. 



2.1 Near-infrared Data 

Observations in the near-Infrared were carried out during the night 
of October 11, 2004 using the Omega 2000 prime focus camera 
(Baumeister et al. 2003) at the Calar Alto 3.5m telescope. Omega 
2000 is equipped with a 2048 x 2048 pixel HAWAII-2 array de- 
tector. It provides a pixel scale of 0.45 arcsec/pixel and a total field 
of view of 15.4 x 15.4 arcmin^ on the sky. Images were obtained 
in the H2 1-0 S(l) line using a 1% narrow band filter centred at A = 
2.122/im, and in a 1% narrow band continuum filter centred at A = 
2.143/im. The total exposure time was 1500s per filter with the av- 
erage seeing of ^ 1.0^^ The data reduction proceeded through stan- 
dard routines using CCDPACK and KAPPA data reduction pack- 
ages developed under the Starlink Project: i.e. the images were dark 
subtracted, flat fielded, sky subtracted and then combined together 
to form a mosaic picture of the complete field. In the next step 
images in A = 2.122/im and 2.143/im wavelengths were astromet- 
rically calibrated and subtracted in order to reveal the pure H2 1-0 
S(l) line emission objects in the field. 

2.2 Mid-infrared Data 

We acquired IRAC (Fazio et al. 2004) and MIPS (Rieke et al. 2004) 
data from the Spitzer Science Archive covering the area around 
GM2-4. The basic flux calibrated images of the Spitzer Science 
Center (SSC) pipeline were used. Cosmetic corrections, astromet- 
ric refinement and the final mosaics were constructed using the 



3 RESULTS AND DISCUSSION 
3.1 Definitions and Nomenclature 

Here we give a brief reasoning on the used nomenclature for the 
protostellar objects referred in the following sections. We follow 
the Gutermuth et al. (2009) source numbering with a slight dif- 
ference aiming at clarifying the overall situation in the immediate 
vicinity of the GM2-4 and IRAS 05373-^2349 objects. Gutermuth 
et al. (2009) call these objects the GGD4 cluster of sources, but 
this nomenclature in our view is misleading in two ways. Firstly 
the actual position of the GGD 4 object reported originally by Gyul- 
budaghian et al. (1978) is about 1.5 arc-minutes to the north-east of 
the near-infrared cluster discovered by Hodapp (1994) (Fig. 1(b)). 
Secondly the GGD 4 object, which were later renamed GGD 4a and 
GGD 4b (GM2-4) by Carballo et al. (1988), are a pair of reflec- 
tion nebulae. These reflection nebulae are not associated with HH 
objects; nor are they obviously connected with the maser sources 
reported by Rodriguez et al. (1980). We believe that the cluster 
should be called the IRAS 05373+2349 cluster and the numbering 
of the infrared sources reported by Gutermuth et al. (2009) should 
be preserved. Therefore in the interest of brevity, and because of 
the absence of any other clusters in the region we will use the IR 

^ GAIA is a derivative of the Skycat catalogue and image display tool, 
developed as part of the VLT project at ESO. Skycat and GAIA are free 
software under the terms of the GNU copyright 
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Figure 3. Area containing the MHO 735, 736 and 737. Positions of IRS ob- 
jects are marked with the crosses on all panels. Dotted arrow on both panels 
show the possible outflow direction. Panel (a) shows the H2 1-0 S(l) emis- 
sion. Panel (b) is a colour composite image constructed from 3.6/im (blue), 
4.5/im (green) and 8.0/im (red) Spitzer IRAC bands where the position of 
MHO 735 is marked with a filled circle for an easy comparison and a dotted 
bow-shock outlines the MHO 737. 

cluster and the IRS# (Infra-Red-Source) nomenclature throughout 
this work. 



3.2 New Molecular Hydrogen emission line Objects (MHOs) 

Figure 2 presents a part of our region where we detected numerous 
H2 1-0 S(l) line emission objects. In order to identify pure H2 1-0 
S(l) line emission objects we subtracted the 2.143/im narrow-band 
continuum image from the 2.121/im line image. For the compar- 
ison we include optical [SII] line data (Nikogossian et al. 2009) 
along with the colour-composite figure constructed from Spitzer 
IRAC bands. We detected number of new Molecular Hydrogen 
emission-line Objects (Davis et al. 2010) (MHOs^ hereafter). Ta- 
bid provides the coordinates, flux in units of 10~^^ W m~^ 
arcsec"^ and some additional notes related to the newly discov- 
ered MHOs. In the following subsections we briefly describe each 
outflow. 



Table 1. List of MHOs detected in H2 1-0 S(l) with cross-identifications 
with HH objects from Nikogossian et al. (2009) and the knots identified by 
Varricatt et al. (2010) are marked by [VDR]8-# nomenclature. 
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t Flux in 10 ^^Wm ^ units and the background Icr noise estimate is 
1.15 X 10~^^ W m~^ calculated in 6" circular aperture 

due to the large point spread function of IRS 3 (cf. Fig. 3b). In con- 



trast MHO 736 and 737 are detected at 4.5/im as expected if the 
emission originates from shock excited H2 lines (Smith & Rosen 
2005; Ybarra & Lada 2009). A confirmation of the shocked nature 
in 4.5/im would require a spectroscopic follow-up. 

A close examination of the colour-composite IRAC view 
(Fig. 3b) reveals the bow-shock shape of the MHO 737 which sug- 
gests the IRS 6 as a possible driving source of a flow with a length 
of at least 0.195pc (Tab. 2). The direction of the flow (if confirmed) 
would also be suitable for inclusion of MHO 736 as a forward 
shock making it over 0.325pc long. 



32. 1 MHO 735, 736 and 737 

About 2 arc-minutes west from GM 2-4 (cf. Fig. 2a,b) three emis- 
sion line objects MHO 735, 736 and 737 were detected. The de- 
tailed view of the area in near-infrared H2 1-0 S(l) emission and 
the composite constructed from the Spitzer IRAC bands including 
these objects are presented on Fig. 3a and b. 

MHO 735 is situated very close to the continuum source IRS3 
which appears prominently at mid-infrared wavelengths and is 
most likely the driving source of MHO 735 based on the positional 
closeness (^ O.OlSpc on 1.17kpc). We were unable to identify 
MHO 735 in the 3.6/im, 4.5/im and 8.0/im Spitzer IRAC bands, 

^ http: //www. jach.hawaii . edu/UKIRT/MHCat / 



3.2.2 MHO 734, 739, 740, 744 and 745 

A detailed view of the region around IRAS 05373-h2349 (IRS4 and 
5) shows a complex structure of knots and bows in the immedi- 
ate vicinity of IRS4 (cf. Fig. 4a). This region has been studied by 
Varricatt et al. (2010) who also detected several faint H2 1-0 S(l) 
line features (object 8 in their paper) which coincide with some 
of the individual knots detected in our study (cf. Tab. 1 for cross- 
identification). In our deeper images we were able to detect more 
features and therefore suggest several possible outflows, which are 
marked on Fig. 4b and listed in Table 2. 

An elongated object MHO 739 located about 10^^ south-west 
from IRS4 was also detected by Varricatt et al. (2010) who la- 
bel this feature as 8-6. This object could be connected with the 
MH0744B compact knot (8-3) and MH0744A bow-shock (8-2) 
to form a 0.144pc long bipolar outflow (739-744 outflow) possi- 
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Figure 4. Close-up view in H2 1-0 S(l) emission of the area surrounding the IRAS 05373-1-2349 and the near-infrared cluster where numerous MHOs where 
identified. Positions of IRS4 and 5 are marked with crosses on all panels. Panel (a) reveals details of MHO 734, 739, 740, 744 and 745. Panel (b) shows the 
suggested outflow structure with dotted arrows. Overlaid contours are from the 3.6cm emission starting from 3cr and climbing with a 3cr step-size and are from 
the data published by MoHnari et al. (2002). Positions of 3.6cm peaks are marked with CM-A, -B and -C. 



bly driven by IRS4/CM-A (Fig. 4b). Next bipolar flow (745flow 
^ 0.068pc), perpendicular to the 739-744 outflow, connecting 
MHO 745A and B on one side and MHO 745C on the other side is 
most likely driven from the very same IRS4/CM-A source (Moli- 
nari et al. 2002). 

The chain of bow-shocks, compact knots and elongated fea- 
tures denoted as MH0734A to G (cf. Fig. 4a) connect with the 
MHO 740B and 740A knots forming the third bipolar outflow (734- 
740 ^ 0.226pc; Tab. 2) in this active region (Fig. 4b). Determina- 
tion of the true driving source of this bipolar outflow is unrealistic 
due to the limiting dataset, but one could suggest some candidates 
based on a simple geometrical placement of the probable candi- 
dates. In this respect CM-B (Molinari et al. 2002) is the most prob- 
able driving source for this outflow but one could also consider 
IRS4/CM-A source because of its apparent closeness and promi- 
nence. Lastly the overall outflow direction described by 734-740 
(Tab 2) coincides with the ^^CO outflow reported previously by Ca- 
soli et al. (1986). 



3.23 MHO 738 and 741 

Figure 5 gives the details of MHO 738 and 741 in H2 1-0 S(l) emis- 
sion and in Spitzer IRAC composite view where the offsets are cal- 
culated from IRS4 (cf. Fig. 2). 

About 40^^ south-west from the IRS4 several bow-shaped fea- 
tures labelled MHO 73 8 A to C are located (Fig. 5c and d). The mor- 
phology and their placement suggests that one of the sources in 
the IR cluster drives an outflow creating the MHO 738 A-C bow 
shocks. To the north-east of the IR cluster, again at a distance of 
about 40'' the two faint knots MHO 741 A and 741B are located 
(Fig. 5 a and b). 

MHO 738 and 741 quite likely are part of yet another bipo- 
lar outflow 0.498pc; Tab. 2) driven from the sources in the IR 
cluster. As a candidate driving source for the 738-741 outflow both 
IRS4/CM-A and CM-B can be considered. 
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Figure 5. Scaled-up views of MHO 738 and 741 where the offsets are cal- 
culated from IRS4. Dotted-line arrows show the direction of the suggested 
bipolar outflow. Panel (a) shows the area including MHO 741 in H2 1-0 
S(l) emission and the panel (b) is the same area in Spitzer IRAC compos- 
ite view constructed from 3.6/im (blue), 4.5/xm (green) and 8.0/im (red) 
bands. Panel (c) and panel (d) show the same wavelengths as in panels (a) 
and (b) but for the area including the MHO 73 8 A, B and C. 



3.2.4 MHO 742 and 743 

The area around GM2-4 has been studied recently by Nikogos- 
sian et al. (2009) who found two Herbig-Haro objects, HH 941 and 
HH942. Our follow-up on the same region reveal near- to mid- 
infrared counterparts of these HH objects. Fig. 6 gives the optical 
[SII] and I-band, near-infrared H2 1-0 S(l) and Spitzer IRAC views 
of the region around the optically undetectable IRS7 which illumi- 
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Figure 6. Region around IRS7 including MHO 742 and 743. Panel (a) is 
optical [SII] view where HH941 and 942 are marked. Panel (b) is optical 
I-band continuum view where the reflection nebulae GM 2-4 (GGD 4b) and 
GGD 4a are indicated. Panel (c) is the near-infrared H2 1-0 S(l) line image 
where MHO 742 and 743 are indicated. Dotted arrow shows the suggested 
bipolar outflow. Panel (d) is Spitzer IRAC composite view constructed 
from 3.6/im (blue), 4.5/im (green) and 8.0/im (red) bands where along with 
MHO 742 and 743 IRS7 is also indicated. 



Table 2. Outflows and their driving sources for the discussed region. 



MHO outflow 


Source 


Size (pc)^ 


Flux* 


Notes 


735flow 


IRS3 


0.018 


13.9 


monopolar 


736flow 


IRS6 


0.325 


15.9 


monopolar 


738-741 


IRS4/CM-A 


0.498 


28.3 


bipolar 


739-744 


IRS4/CM-A 


0.144 


63.1 


bipolar 


745flow 


IRS4/CM-A 


0.068 


41.3 


bipolar 


734-740 


CM-B 


0.226 


165.3 


bipolar 


742-743 


IRS7 


0.372 


104.1 


bipolar 



t Assuming 1.17kpc distance 

~^ Combined flux in 10~^^ Wm~^ units 



nates the GM 2-4 (GGD 4b) nebula. For completeness we note that 
GGD 4a is also illuminated by a near- to mid-infrared source which 
was classified as a Class II pre-main- sequence star by Gutermuth 
et al. (2009) (cf. Fig 1). 

About 20'' south-west of IRS7 (Fig. 6c) the MHO 742 object 
composed of several prominent knots (A, B and C) was identi- 
fied coinciding with the optical HH 941A and B (cf. Tab. 1). While 
knots MH0742B and C are less defined the MH0742A knot 
clearly resembles a bow-shock like structure pointing back to IRS7. 
Similarly about 30'' north-east of IRS7 (Fig. 6c) lays MHO 743 
which has a clear bow-like morphology pointing back to the IRS 7 
as well. This makes 742-743 a 0.372pc long bipolar outflow (Tab 2) 
emanating from IRS 7. 

In addition all features detected in near-infrared H2 1-0 S(l) 
line are also clearly detectable in green on the Spitzer IRAC com- 
posite view (Fig. 6d) which suggests that the radiation must be 
coming from the shock excited H2 lines found in IRAC 4.5/im 
band. 



3.3 SED modelling of IRS Objects 

In order to better characterise the IRS objects we have explored 
the archive of 2D axisymmetric radiative transfer models of proto- 
stars (Robitaille et al. 2006) calculated for a large range of stellar 
parameters (masses, accretion rates, disk masses and disk orienta- 
tions), evolutionary stages (from the early envelope infall to the late 
disk-only) and each calculated for 10 different viewing angles de- 
veloped by Robitaille et al. (2007). The archive also provides an 
SED fitting interface^ which can be used to select best fit models 
to a given observed SED within the specified ■ We have used 
this interface to input all the available magnitude/flux values for 
the IRS objects from various archives and databases (presented in 
Table Al). Those values were imputed assuming a conservative er- 
ror of at least 10% in order to avoid over-interpretation of the SEDs 
(Robitaille et al. 2007), although the true errors are smaller in the 
near to mid-infrared range. The observed SED can also be scaled 
to match the models in the grid by using a range of values in ex- 
tinction (Ay) and distance. For the distance we took range from 
1 to 1.7 kpc but for IRS4 and 5 we narrowed down the range to 
1-1.3 kpc, since the distance to the IR cluster (1.17 kpc) is known 
with higher certainty. For the Ay we used values from Gutermuth 
et al. (2009) with zb 10% range. 

Figure 7 shows the SED fits for the IRS objects where the solid 
black line in each case represents the best fit model and the grey 
lines show the next best fits which satisfy the criteria of - xLst 
< 3 per data point (Robitaille et al. 2007). In the case of IRS4 we 
present a second set of models which better fit the long wavelengths 
part of the SED (IRS4 long). Results from the fits are then listed in 
Table 3 where first column lists the names of fitted parameters and 
in all subsequent columns the parameters of SED fits are given for 
each IRS object. 

Before discussing each object and determining how accurate 
the derived physical parameters are from the given SED fit, we 
would like to address some of the general discrepancies that one 
could expect from the procedure itself. We emphasise that the SED 
fitter is not a tool to accurately derive the physical parameters but 
it provides a way to determine how well constrained each of the 
different parameters are (Robitaille et al. 2007). Those parameters 
strongly depend on the wavelength coverage for each object and 
are better constrained by the multi- waveband data points from the 
near-infrared to radio (see Table Al for our coverage). 

In cases when the observed data-points were enough to fit a 
relatively accurate SED there have been some cases (CoKuTau/1 
and GG Tau for example) when the derived parameters where quite 
unrealistic when compared with the observed data and the indepen- 
dent theoretic analysis (see Robitaille et al. 2007, for the details). 
The SED fitter uses pre-compiled grid of model SED (Robitaille 
et al. 2006) to derive the physical parameters and as numerous as 
they are it would be naive to suggest that we have a complete cov- 
erage of all the possible nuances that influence or drive the star 
formation process. We would like to note that it is not the scope of 
this work to give exhaustive explanation to all those discrepancies 
but rather to use these SED fits as first approximation or as a work 
in progress. 



^ http: //caravan . astro .wise . edu/protostars/ 




Figure 7. Spectral energy distribution (SED) of IRS 3 to IRS7. Filled circles are the data points and the triangles are lower and upper limits where errors are 
generally smaller than the data points. The solid black line indicates the best-fitting model, and the grey lines show all models that also fit the data satisfying 
the - xlest ^ ^ ^^^^ point criteria. The dashed line shows the SED of the assumed stellar photosphere in the best-fitting model. 



Table 3. Model best fit results (in bold face text) satisfying the - xlest ^ ^ ^^^^ point criteria for IRS objects using online SED fitting 
interface by Robitaille et al. (2007). Provided errors are calculated as difference between lowest and highest parameter values of SED models 
compared to the best fit satisfying the same criteria. See Section 3.3 for explanation and discussions. 



Parameter 


IRS3 


IRS4 


IRS4-long 


IRS5 


IRS6 


IRS7 


M* (Mq) 




9 90+2.10 

^•^"^4.64 




3 55+4.78 

•^•""^3.16 


0.29±3:« 


2 24+0.00 
^•^^^0.00 


Age X 10^ (yr) 


4 83+0.00 


138±250 


0.58±J:f, 


"•^^^0.10 




54 6±o.o 


R. (R©) 


10 6±0-0 






25 7+35.0 


4 75+2.54 

^•^"^3.01 


4 93+0.0 


T* X 10^ (K) 


43 1 0.0 


25.5±iJ 


^•^•^^0.12 


4 21+15.09 
^•^^^0.52 


3 35+9.85 


4.96±g:g 


Mdisk X 10-2 (M0) 




7.55±ioj5 




20.40±i9.io 






Md^sk X 10-6 (M0/yr) 


0.024±0:00 


0.067±6.333 






0.081±97.9i 




Rmin Disk (AU) 








5.89^1^21^ 


57+2.32 
"•"^^0.62 


0.33±°:° 


Rmax Disk (AU) 


34.8±°:° 


285±208 




11 q+553.1 


AQ CI 2951.5 
^0-^^46.75 


83.0±°:° 


Menv. X 10-2 (Mq) 


5 57^0.00 

"•''^^0.00 


0.0000285±«:g099p 


572+13728 


100±2|80 


35.8±35;g 


"•"'=■=0.00 


Menv. X 10-6 (M0/yr) 


o iQ^O.OO 
^•-^^±0.00 




454+1296 




7.5±?^5^ 


"•"'=■=0.000 


Rmax Env. xl03 (AU) 






5.04±69^56 


8.7±|9,^ 






Ltot (Lq) 


34.5±o.o 


/:AaA+5570 
""•^"^5450 


290±i96 


212±?ei 


2.68±i27^32 




Incl. Angle (deg.) 


56.6±6.7 


81-4±!o'o 


18.2±o.o 


31.8±«:| 


49.5±37.6 


31.8±24^8 


Fitted Ay (mag.) 


22.00±5;00 


22.68±3.32 


22.00±J:56 


38.76±ii:i| 


16.3±3.7 


14 91+0.00 


Fitted Distance (kpc) 


1 15+0.00 


1 12+0.17 


1 21+0-06 


1 17+0.11 

'•"='=0.17 


1 55+0.15 


1 17+0.05 

-*^--*^^^o.oo 


N (models) 


2 


19 


24 


78 


96 


4 



IRS 3 was previously reported to be a Class I protostar by Guter- 
muth et al. (2009) using near to mid-infrared colour-colour analy- 
sis. Our SED modelling gives a good fit (Fig. 7) which is largely due 
to the numerous multi- waveband data points (1.2 to 70 /xm) present 
for the object. The degeneracy or number of models satisfying the 
best- fit criteria is just 2 with slight difference in disk inclination 
angle, fitted distance and the fitted Ay (Tab. 3). The fitted distance 



(^l.lSkpc) is in good agreement with the value inferred from the 
observations of the IR duster (Molinari et al. 1996). The same 
can be said about the fitted Ay (^^20-22) which is also in a good 
agreement with the value derived by Gutermuth et al. (2009). The 
error in fitted inclination angle is not unexpected since this is one 
of the least well determined parameters in the SED fitting process 
but it is directly connected with the discrepancy in distance and Ay 
determination which is the result of an input ±10% assumed range. 
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The rest of the fitted parameters do not show any degree of 
degeneracy. However, due to the fact that our data does not ex- 
tend beyond the Spitzer MIPS TO/xm we suspect that some of the 
parameters are not well constrained. In particular the absence of 
any (sub)mm data points makes the precise estimates of the disk 
mass {Mdisk) questionable and introduces a degeneracy between 
the central source temperature (T*) and the disk accretion rate 
(Mdisk)- On the other hand we can stipulate that SED fit can not be 
interpreted by disk-only source but there is a significant rise after 
^lOfim which can only be interpreted by the presence of an enve- 
lope and its higher rate of accretion (Menv.) (Robitaille et al. 2006). 
This last point strongly supports the protostellar nature of the IRS 3 
object confirming the findings by Gutermuth et al. (2009). The clas- 
sification of this source and the positional placement in regards to 
MHO 735 leaves no doubt that IRS3 is the driving source of a short 
(0.028pc) outflow. 

3.3.2 IRS4 and IRS4 long 

IRS4 is one of the prominent members of the IR Cluster. It has 
been observed at multiple wavelengths between 1.2 to 1300/im (as 
can be seen from the fitted SED (Fig. 7). The prominence of this 
object enabled its early discovery as a red source called CPM 19 
which was then resolved into IR Cluster (Campbell et al. 1989; 
Carballo et al. 1988; Hodapp 1994), as noted in Section 1. Recently 
Nikogossian et al. (2009) observed the same object in optical I and 
R bands and reported about the variability comparing DSSl, DSS2 
and current data. In addition to these facts IRS4 was classified as a 
Class I protostar by Gutermuth et al. (2009) based on near to mid- 
infrared colour-colour analysis. 

The observational facts would seem to suggest that IRS4 is 
quite an energetic source and it would be misleading to place it into 
one or other classification group. Long- wave observations suggest a 
possibility of high-mass source (Molinari et al. 2008) which would 
make low-mass classification of YSOs incompatible or inaccurate. 
It is worth mentioning that all of the previous infrared to radio ob- 
servations were done using larger apertures which would invali- 
date any kind of comparison with the current near to mid-infrared 
data due to the IR Cluster source contamination. With these facts 
in hand we performed two separate SED fittings of IRS4 object; 
in one SED fit we used only near to mid-infrared data obtained 
with a small aperture (^5000AU) and the other we used large- 
aperture (^22000AU) data to incorporate the flux values (such as 
MSX6, c.f. Tab. Al) from lower resolution mid-infrared to radio 
wavelengths (IRS4 long). Taking into consideration all the possi- 
ble discrepancies in the data points it would be almost impossible 
to correctly derive any physical parameters. However, as a first ap- 
proximation we tried to give some interpretations noting the highly 
speculative nature of them. 

When we fit an SED to the near to mid-infrared data points 
(Fig. 7) we get edge-on system containing zero age main sequence 
star or ZAMS of about 1OM0 with almost no envelope and rela- 
tively small accreting disk. The fitted temperature would seem to 
indicate a spectral type of late O or early B with total luminosity 
of 6000 L0 which would explain its observability in all the wave- 
lengths. In contrast, if we fit the data-points represented by larger 
aperture observations (Fig. 7) we get about 4.5M0 Class O/I ac- 
tively accreting protostar with a large envelope which would seem 
to fit the role of outflow-driving source. Both of these fits have con- 
siderable level of degeneracy which is understandable due to the 
uncertainties in the apertured used (Tab. 3). 

This kind of clear separation between sources fitted using the 



different apertures suggests a possibility of multiple objects in the 
line of sight belonging to the same IR Cluster. It would seem that 
the visible IRS4 is situated in front of an embedded core which 
probably contains Class O/I protostar coinciding with the CM-A 
core (Fig. 4). This kind of superimposition would explain the abun- 
dance of outflows from the same area but it would be difficult to 
stipulate which outflow belongs to which source since there is a 
possibility of more hidden objects along the line of sight. 

3.3.3 IRS5 

IRS5 is one of the Class I sources identified by Gutermuth et al. 
(2009) which is almost invisible at near-infrared wavelengths. The 
picture is quite different in mid to far-infrared bands where IRS5 
even outshines IRS4 (24/im on Fig. 1 ; Tab. Al). Our analysis of this 
object provides a challenge due to its apparent closeness to IRS4, 
preventing accurate PSF flux measurements. Similar discrepancies 
are also invalidating some of the flux values presented in various 
catalogues forcing us to use them only as upper-limits which rises 
the degeneracy of the SED fits (Tab. 3). Despite that, from the SED 
fit (Fig. 7) we can derive some reasonable characterisation for IRS5. 

IRS5 is one of the deeply embedded members of this IR Clus- 
ter of stars. Our SED fit suggests a very young Class O/I protostar 
possibly in the same stage as CM-A core discussed in Sect. 3.3.2. 
There was no 3.6cm core reported in this position (cf. Fig. 4) within 
3a limit of the observations (Molinari et al. 2002) but there is a hint 
of a possible object below or around that limit which would need 
to be followed by separate observations. IRS5 could be the driving 
source of the MHO 740 A and B flow but those MHOs seem to be 
aligned with the knots in MHO 734 and are likely being driven from 
the CM-B core rather than from IRS 5, which is not well aligned 
with the axis of this MHO outflow (Fig. 4). 

3.3.4 IRS6 

One of the possible candidates of outflow driving sources in the pre- 
sented region is IRS6. This object was identified as a Class I YSO 
by Gutermuth et al. (2009), based solely on its mid-IR colours (it 
is undetected in the near-IR). The scarceness of valid data points 
throughout the whole electromagnetic spectrum prevented us from 
producing a reliable SED fit for IRS6 as can be seen from Fig. 7 and 
Tab. 3. The best fit SED values seem to indicate a young Class O/I 
low-mass star (O.29M0) situated relatively further (1.55kpc) than 
the rest of the IRS objects being discussed here. All of the fitted 
parameters in this case are very sensitive to the range of Ay used 
which is still uncertain for this particular source. In our fit we as- 
sumed range of 14 - 20 based on the fact that values of 14 to 20 
were reported for IRS 7 and IRS 3 respectively by Gutermuth et al. 
(2009) and IRS6 is situated in-between. 

3.3.5 IRS7 

IRS7 is the star which illuminates the GM 2-4 (GGD 4b) optical re- 
flection nebula and is the probable driving source of the MHO 742- 
743 outflow (Fig. 6). We have relatively well sampled data points 
from near to mid-infrared range of electromagnetic spectrum which 
enabled us to produce quite good SED fit as can be seen on Fig. 7. 
Following the same logic as we did in case of IRS 3 (see Sect. 3.3.1) 
we can suggest that the best fit SED parameters (cf. Tab. 3) infer 
relatively evolved Class II YSO rather than Class I as suggested by 
Gutermuth et al. (2009). We find that IRS7 is a low-mass YSO 
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Figure 8. The wide-field view of the region in H2 where offsets are mea- 
sured from IRS4 and solid line box shows the area presented in Fig. 2(b). 
Overlaid contours are from map constructed from 2MASS PSC and 
provided by Rowles & Froebrich (2009). Filled circles indicate the posi- 
tions of IRAC identified Class I/II YSO sources provided by Gutermuth 
et al. (2009) where the colour spans according to the range of IRAC SED 
Slope (a^^^^). 

with active accreting disk which is dominating the emission com- 
ing from this source. All the fitted SED parameters enforce our 
scenario of IRS7 being the driving source of the 742-743 bipolar 
outflow. 



3.4 Star Formation in the vicinity of GM 2-4 

Fig. 8 shows the wide-field area where the search for outflows was 
performed. However, MHO flows were only detected in a small 
sub-region presented in Fig. 2(b). This suggests that the immedi- 
ate vicinity of the IRS3 to IRS7 is the site of relatively recent star 
formation activity. This scenario becomes even more favourable if 
we consider the disk prominence of the YSOs as defined by the 
parameter a^^"^^ (the slope of the SED across the mid-infrared 
wavelengths observed with the IRAC instrument on Spitzer) and 
their overall distribution in East and West of the area. There are 
more Class II or weak-disk sources in the East than in the west. 

Furthermore let us consider a^^^^ and the relative distribu- 
tion of the YSOs to the two distinct Ay peaks (A and B) shown 
on Fig. 8. Most of the YSOs with significant circumstellar disks are 
associated with the Ay peaks but not strictly confined to the peaks 
which could generally be due to the dispersal of protostars from 
their original birthplace showing somewhat uniform spread in rela- 
tion to the A^ peaks as it can be seen in case of peak A. This would 
mean that YSO population associated with peak A could have al- 
ready gone through the active phase of star formation and therefore 
had time to disperse almost uniformly in relation to the peak A and 
have more evolved appearance. 

The same could not be said for the YSO population associ- 
ated with the peak B (Fig. 8) since in this case there is no uniform 
distribution rather than the protostars are situated in the North-East 
part of the peak suggesting that the star formation progresses from 



North-East to West, South-West for the entire region discussed 
here. This is in good agreement with the scenario of star formation 
around Ay peak A and non-detection of outflows in that area. 



4 SUMMARY 

We have presented near-infrared H2 1-0 S(l) line study of the re- 
gion surrounding GM 2-4 alternatively known as GGD 4 in order to 
reveal the outflow content of the area and to better understand the 
earlier suspected star formation processes in this area. In our quest 
to understand this region we combined multi- waveband data- sets 
obtained from various online archives. Our findings and analysis 
can be summarised as following: 

• The analysis of our near-infrared H2 1-0 S(l) line data enabled 
us to identify 27 individual shocked features which combine into 5 
bipolar and 2 monopolar outflows driven from the protostars in the 
immediate vicinity of the GM 2-4 reflection nebula. We were able 
to cross-correlate some of the H2 identified objects with previously 
reported optical HH objects as well as with green fuzzy features 
in Spitzer IRAC 4.5/im wide-field band thought to have a shocked 
nature. 

• We were able to propose the driving sources for the detected 
outflows based on the morphological and positional characteris- 
tics of the individual shocked features. All of the driving sources 
were identified as YSO based on the previously available knowl- 
edge (Gutermuth et al. 2009; Nikogossian et al. 2009) and our SED 
modelling. In one particular case, IRS4, our SED modelling en- 
abled us to suggest a solution for the previously known discrep- 
ancy (Molinari et al. 1996, 2002) as two different stars one late O 
to early B type 1OM0 ZAMS high mass star and the other as very 
early, possibly Class protostar on the same line of sight which 
might be the actual driving source of the several bipolar outflows. 

• Our analysis of the wide area using IRAC SED Slope in con- 
junction with the A^ map including the GM2-4 suggests that 
star formation progresses roughly from North-East to West, South- 
West. This result confirms the conclusion made in our previous pa- 
per (Nikogossian et al. 2009) mainly on the base of the optical data. 

There are still several facts that we wish to follow-up in or- 
der to understand the exact outflow dynamics involved in the IR 
cluster for which further high-resolution imaging, spectroscopy 
and radio-line data would be beneficiary. The importance of spec- 
troscopic followup and monitoring of IRS4 (CPM19) would prove 
useful in determining the exact nature of variability in this object. 
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APPENDIX A: PHOTOMETRY 

This paper has been typeset from a TgX/ LT^X file prepared by the 
author. 
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Table Al. Photometric value database for all IRS object converted into mJy with the appropriate errors and used apertures to model SEDs. Origin of the each 
data-set is indicated in the last column. 



Filter 
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^OS 747-1-^ fx^fx (A\ 


^ 01 ^-1-0 0^^^ (^A\ 
j.y i-jzizyj.yjjKj 


f\d 0/^7-1-0 ^Qf\ 




TP AP PVizl 




90^1 ^^0^-1-9^ 874 ('^\ 


Id^ 194-1-^ 81^ fA\ 


4 ^^»1 -1-0 1 94 fA\ 


81 4S9-I-0 747 ('^\ 








^91^»0 7 r^l S'i 








(A\ 


Al^APT/TPP Q 


89-1-0 99^ ('^\ 


^tn.^ozizi /J.^ \o) 








\?) 






^0090 ^"^0^ 








yp) 


MSX6 P 




<r46S2 9 H 














^70QS ^ n S^i 
^ i\jy Omj v^J^ 








(A^ 


AKART/TRP 1 R 


COO c^icy g (f\ 


TT499 8+5^ 4 








\-^) 


iVlwjyvO Cj 




^19474 n 








(A\ 


MIPS 24 


1174.08±31.99 (7) 


> 6768.59 (4) 


> 6945.42 (4) 


33.55±1.21 (5) 


243.05±2.23 (7) 


(2),(3) 


IRAS 25 




<305560 (30) 








(6) 


IRAS 60 




< 144830 (30) 








(6) 


AKARI/FIS 65 




115140±7332 (30) 








(7) 


MIPS 70 


<3223 (30) 


<96926 (30) 








(3) 


AKARI/FIS 90 




< 120067 (30) 








(7) 


IRAS 100 




<191160 (38) 








(6) 


AKARI/FIS 140 




132767±23968 (30) 








(7) 


AKARI/FIS 160 




239313±3880.5 (30) 








(7) 


350/im 




<41000 (30) 








(8) 


SCUBA 450/im 




< 18000 (30) 








(8) 


SCUBA 850/im 




<2730 (35) 








(8) 


1.1mm 




<910 (30) 








(8) 


1.3mm 




<700 (30) 








(8) 



t Data origin: (1) 2MASS by Skrutskie et al. (2006); (2) Gutermuth et al. (2009); (3) this publication; (4) MSX6 PSC by Egan et al. (2003); (5) AKARI/IRC 
by Ishihara et al. (2010); (6) IRAS PSC; (7) AKARI/FIS by Kawada et al. (2007); (8) Molinari et al. (2000). 



